Euclidean 4d exact solitons in a Skyrme type model 
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Abstract 

We introduce a Skyrme type, four dimensional Euclidean field theory made of 
a triplet of scalar fields n, taking values on the sphere 5^ , and an additional real 
scalar field <j), which is dynamical only on a three dimensional surface embedded 
in ]R . Using a special ansatz we reduce the Ad non-linear equations of motion 
into linear ordinary differential equations, which lead to the construction of an 
infinite number of exact soliton solutions with vanishing Euclidean action. The 
theory possesses a mass scale which fixes the size of the solitons in way which 
differs from Derrick's scaling arguments. The model may be relevant to the study 
of the low energy limit of pure SU{2) Yang-Mills theory. 



The aim of this Letter is to present a new four dimensional Euchdean field theory 
and to construct an infinite number of exact soliton solutions for it. The physical moti- 
vation of our work is twofold. First, it contributes to the development of exact methods, 
in the framework of integrable field theories, for the study of non-perturbative aspects 
of physical theories [T]. Second, it may connect to recent attempts to understand the 
role of solitons in the low energy limit of pure Yang- Mills theories |2]. 

The Skyrme Pj and Skyrme-Faddeev jlj models are effective field theories possessing 
topological solitons, with several applications in various areas of physics. Such theo- 
ries can not be solved exactly, and those solutions are constructed through numerical 
methods. The lack of an exact closed form for the solitons prevents the understanding 
of several properties of the phenomena involved. It is therefore of physical interest to 
construct similar theories possessing exact solutions. In Ref. [Sj we have introduced 
a model with the same field content and topology as the Skyrme-Faddeev theory, and 
have constructed an infinite number of exact soliton solutions with arbitrary topologi- 
cal Hopf charges. That model is integrable in the sense of [Ij and possesses an infinite 
number of local conservation laws. The exact solution was obtained through an ansatz 
that, as explained in jH|, originates from the conformal symmetry of the static equations 
of motion. 

In this Letter we propose a model with a Skyrme-Faddeev scalar field taking values 
on the sphere 5^, together with an extra real scalar field. An important aspect of the 
model is that, although defined on a four dimensional Euclidean space IR'^, the extra 
field lives only on a three dimensional surface embedded on IR^. That fact introduces 
a mass parameter which fixes the size of the solitons, but in a way which differs from 
Derrick's scaling arguments. The soliton solutions are constructed through an ansatz, 
introduced in [UllH], which reduces the four dimensional non-linear equations of motion 
into linear ordinary differential equations. The physical boundary conditions are such 
that the Euclidean action vanishes when evaluated on all soliton solutions. The theory 
can perhaps be viewed as an effective theory for the pure SU{2) Yang-Mills theory in 
a way explained at the end of the Letter. 

The model considered in this Letter is defined by the action 

h (1+ I u |2)2 e2 ilR^ f"" ^ ' 

where and u are real and complex scalar fields respectively, with the latter related 
to the triplet of fields n living on S"^ {n^ = 1) through the stereographic projection 



n 



(u + u, -I {u-u),\u |2 -l) / (l+ I u 1^) (2) 



H^y is the pull-back of the area form on S"^ 

^ _^ j^a,u3u- - 3ua,u-) ^ . ^ g^ _j J3J 

(1+ I U Y) 

M and e^ are coupling constants, with M having dimension of mass and e^ being di- 
mensionless. E is a 3-dimensional surface embedded in the four dimensional Euclidean 
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space-time H^. As we explain below E corresponds to the surface where a special 
variable y vanishes, and so one can write /x; c^S = j^id!^xuj 5 {y)^ with uo being an 
integration measure defined below (see (jT^). 

The action (^Q) is not invariant under the 5*0(3) transformations on the sphere S'^ 
defined by the fields n. The first term of (^ breaks that symmetry. However, (^ is 
invariant under the global U{1) ®U{1) defined by the transformations 

u -> e*" M 0^0 + /? (4) 

The Euler-Lagrange equations associated to (dj) are 

d^ {H.^d'^u) +'-Me'u6 {y) ^^^^ = (5) 

together with its complex conjugate, and 

V(i+ I « r) / 

Following [3 El we introduce the ansatz 



based on a set of orthogonal curvilinear coordinates {y, z, ipi, ip2), with — oo < y < oo, 
0<z<l, 0<{pi< 2tt, i = 1,2. We have that g = g {y, z), < g < 1, with m,, 
i = 1,2, being arbitrary integers, and Ei = ±1, arbitrary signs. 

There are in fact two actions of the type (0) with different surfaces S, and as 
consequence two sets of coordinates {y, z, ipi, ip2), admiting exact solutions through the 
ansatz ((Zj). They are 4d generalizations of the 3d toroidal and 2d polar coordinates. 
The toroidal like coordinates are defined as 

^0 r- ro 



xi = — v^cosv92 a^3 = — vT—z sinipi 
p p 

X2 = — yz sm (p2 X4 = — smh y 
p p 



with p = coshy — a/1 — zcosifii, and tq is a constant with dimension of length. The 
polar like coordinates are defined as 



Xi = ro e^ a/z cos V52 xs = ro e^ Vl — z cos ipi 



X2 = tq e^ \fz sin (^2 a^4 = ''^0 ^ VT—z sin ipi (9) 

The metric is given by 

r^ / dz"^ \ 

ds^ = -| idy'' + ^^ + (1 - ^) difi' + zdip2^\ (10) 



with T = p in the toroidal case, and r = e ^ in the polar case. The surface S is defined 
by the equation y = 0, and in the toroidal case is 

St = {x^ G M^ I a;4 = or r^ ^ cx)} (11) 

with r"^ = x\ + x\ + x\ + x\. In the polar case we have 

Sp = 5^ = {x^ G IR' I r^ = rl) (12) 

The integration measure in the toroidal and polar cases are respectively 

^T = p'^/ro , ujp = e~'^y/ro (13) 

Due to the orthogonal character of the coordinates and the form of the ansatz (0) 
it turns out that the gradient of (p is orthogonal to the gradients of | u | and y. In the 
toroidal case one can easily check that cud'^cj) = —d^ujd'^''(f) = —1e\ rriip^ sin (fi/ (tq v^l — z 
In the polar case one has ujd'^(j) = —d^ujd'^(f) = 0. Therefore, in both cases the equation 
of motion ® for is automatically satisfied by the ansatz (|7j). 

Replacing the ansatz ((7j) into the equation of motion (0) one gets the same equation 
for the profile function g {y, z) for the two types of coordinates (jHl) and ©• It is a linear 
partial differential equation given by 

4^^^9. {nd,g) + dlg+^^5{y) g = (14) 

with Q = ml z + TTi^ {1 — z). It can be solved by separation of variables 

g{y,z) = k{y) h{z) (15) 

giving the linear ordinary differential equations 

k" {y) - X' k (y) + -p 6 {y) k (y) = (16) 

and 

zil-z) d, {{q^z + (1 - z)) d,h (z)) + ^ {qh + (1 - z)) /i (z) = (17) 

where A^ is the separation of variables constant, and 

I ^7ii I , , 

(i = -r^ (18) 

I ^2 I 

We must have < g < 1, and so the solutions of (jT^ and (jTTj) should not diverge. 
Since any solution of those equations can be rescaled by a constant, we can use any 
finite solution with definite sign, to construct a g satisfying that requirement. 

Using the integral representation for the Dirac delta function one gets that the 
solution for (fT^ is 



k{y) = k{0)--^ du^-- (19) 
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From Cauchy's theorem one gets 



oo p^vy 



g-a TT 



Therefore, the consistency of (J19p and (j2(Jj) at y = fixes the value of A to 

lAh^ (21) 

and the solution for k takes the form 

k{y) = k (0) e-^^"''^«l^l/^'^ (22) 

Notice that both integrands in (^ are positive definite. The sign of the measure of the 
S integration is dictated, according to (fT^ . by the sign of tq. Therefore, the relative 
sign of the two terms in ([Q) is given by the sign of (— Me^ro), which due to (J2H) must 
be negative. So, the Euclidean action ((T)) can not be chosen to be positive definite. 

The equation for h (fTTjl is what is called a Heun equation j2j. Having four regular 
singular points, it is a generalization of the Gauss' hypergeometric equation. Indeed, 
for g = 1 eq. (fTTj) reduces to the hypergeometric equation. In our case, the fourth 
singular point will not be a concern since {q'^z + (1 — z)) can only vanish at z = 1 for 
g = 0, and we will not treat that case. 

The boundary conditions that h (z) has to satisfy will be dictated by the require- 
ment of finite action solutions. Replacing the ansatz (jZj) into the action (^J one gets 
for both coordinates (jHl) and that 

S = -^ /" * /' dr.^IL^ L(l- .) iA,f + i3,,f - ^ S (,) A 

e^ J-oo Jo z{l — z) \ 8 J 

(23) 
Using (|T^ multiplied by g, integrating by parts and using the solution (J22|) . one gets 
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S = -^^^^ H ^ ^' 1-1 -^2 h h' |.=o) (24) 

If a solution of ()17|1 does not vanish at z = or z = 1 then its first derivative has a 
logarithmic divergence there. So, in order to have finite action one needs h to vanish 
at 2; = and z = 1. In such case, the action vanishes. In addition, since we need 
< (7 < 1, and since the solution for k given in (j^^ has a definite sign, it follows that 
h can not change sign between z = and z = 1. Therefore, the boundary conditions 
required for h are 

h{0) = h{l) = 0; h{z)^OioTO<z<l (25) 

Notice that (J17p is invariant under the joint transformations z ^^ {1 — z) and q ^-» 1/g. 
Therefore, if hg (z) is a solution of (jl7|) for a given value of q, then 

h/g (z) = hgil- z) (26) 



is also a solution of ()17p for the inverse value of q. We will then construct solutions for 
(I17p in power series around z = Q for < g < 1. The solutions for g > 1 are obtained 
from those through ^I^ . Replacing h {z) = J2'^=o Cn^"' into (fT7|) we get that Cq = and 
Ci is arbitrary. The remaining coefficients are obtained through the recursion relation 

{q' - l) f ^ - (n - 2)^) c„_2 (27) 



n{n — 1 



_ _ (^ _ 1) (^ _ 2) + [q' - Ij (n - lY\ c„_i 

Therefore, h{z) vanishes at 2; = 0, as required by (J2S1)- For a given value of q the 
locations of the other zeroes of h depend upon the value of A. We then have to tune 
A in order for the next zero to fall exactly at z = 1. The remaining zeroes will then 
happen for 2; > 1, as required by (J2S1)- That gives a discrete spectrum for A and the 
size of the solution is fixed by the length scale tq through ()21|) . For g = 1 the solution 
for such a problem is exact, since the series solution truncates, and we get 

h,=i{z)=Az{l-z) K=i = ^ (28) 

For other values of q the series does not truncate and we can tune A numerically. The 
values of A, for some values of g, such that /i (1) = is given in Table [H Due to (j^EI) 
the values of A satisfying the boundary conditions (j25p are the same for q and 1/q. 

As we commented below (|7j) the profile function g has to take values between zero 
and unity. Therefore, we take fc (0) = 1 in (j22|) . In addition, we normalize h [z) such 
that its maximum value, on the interval < z < 1, is unity. We then write 



h,{z) = c{q) [z + Y,Cnz''\ (29) 




with c (g) given in Tabled The c„'s are determined from (j?r|) with cq = and ci = 1. 
The position Zc of the maximum, i.e. the zero of the derivative {h' {zc) = 0), is also 
given in Tabled The form of the functions (j^^ is exemplified in the Figure ^ for the 
case g = 1/6. Consequently, the solutions for the profile function g become 

g,=, (y, z)=Az{l-z) e-2^l^l g, (y, z) = h, {z) g-l^^'')!!^! (30) 

with A fixed by (PT|) and the boundary conditions (J2HI)- Values of A for some g's are 
given in Tabled 

Since the solutions (jHUj) satisfy the boundary conditions ((23), the action (P), eval- 
uated on them, vanishes (see plj) ). However, both terms of the action are finite and 
cancel each other. We can then introduce an energy function E, which is the same for 
the two coordinates systems (jH)) and 0, and it is given by 

E = M f dE-^^^^^^ = \ f ^^4 i^2 327r2 



Q 


AV4 


Hq) 


Zc 


c{q) 


1 


2 


7.542472 


0.5 


4 


9/10 


1.99778564857563(7) 


7.568858 


0.5131694374 


3.8013381236 


4/5 


1.99015726516381(9) 


7.661008 


0.5278864055 


3.6045699997 


3/4 


1.98376526603091(4) 


7.739752 


0.5359452134 


3.5064351049 


2/3 


1.96834721693012(0) 


7.935672 


0.5506341233 


3.3423893780 


1/2 


1.91364355092039(7) 


8.707347 


0.5862821184 


3.0065092866 


2/5 


1.85959184446300(8) 


9.612200 


0.6134572404 


2.7934885325 


1/3 


1.81146288327809(9) 


10.571742 


0.6351236762 


2.6427573006 


1/4 


1.73364482996707(3) 


12.547803 


0.6678866103 


2.4390777756 


1/5 


1.67496824177613(6) 


14.538513 


0.6917764653 


2.3049609403 


1/6 


1.62948847681927(1) 


16.519869 


0.7101532977 


2.2084486426 


1/10 


1.51764973827612(1) 


24.281837 


0.7557187711 


1.9890631347 


1/100 


1.25301610750(6) 


177.9(3) 


0.8719052146 


1.5135800082 


1/1000 


1.1610118(7) 




0.9168161751 


1.3465808941 



Table 1: Numerical results for the solutions: i) Some chosen values of q (see definition 
in ()18|) ) are given on the first column; ii) On the second column are the values of A 
such that hq vanishes at 2; = 1 (see below (J2Z|)); Hi) A (g), given on the third column, 
is defined in (IS^ : iv) z^. given on the fourth column is such that h'^^Zc) = and 
hq{zc) = 1; v) c{q) given on the fifth column is the normalization of hg (see ()29|) ) such 
that its maximimum value is unity on the interval < ^ < 1. The numbers are exact 
up the digit shown, except for those cases where the first uncertain digit is shown 
between parenthesis. 



with 



A(g) 




A 




(32) 



The values of A (g) for some of values of q are given in Table Notice that E does 
not depend upon the signs of the integers rrii, and it is invariant under the interchange 
mi <-> r?7,2. 

The solitons we have obtained in this Letter are smooth solutions with vanishing 
Euclidean action, localized around the three dimensional surface S. That surface plays 
the role of a source for the solution. The eq. (0) with the delta term, indicates that 
character. In addition, since the 0(3) symmetry is broken to U{1) f^(l), the theory 
([Q) can perhaps be seen as an effective theory. There can be many applications in 
that direction, but perhaps one that is worth pursuing is the following. Consider pure 
SU{2) Yang-Mills theory (without matter) and the Faddeev-Niemi |2] parametrization 
of the physical degrees of freedom of the gauge field (a modification of Cho's ansatz 

Ii) 



A. 



C^n + d^n An + p d^n + a d^n A n 



(33) 



where p and a are real scalar fields, C^ is an abelian gauge field, and n is the same 
as in (0). The arrows stand for the orientation of the fields in the algebra of SU{2). 



One can use the scalar fields to construct a complex field $ = p + icr =| $ | e*''^. It 
could happen that as the energy decreases the degrees of freedom associated to C^ and 
I $ I get frozen. The phase of $ however gets frozen at different values in different 
regions of space. On the border of such regions it could still have some dynamics. We 
would then have at low energies, the degrees of freedom of the Yang-Mills field being 
described by the field n, all over the space, plus the phase on the border of those 
regions. Very probably the effective action describing such regime is very complex, but 
the action ^ could perhaps describe, in Euclidean space, some aspects of such phase 
of the Yang-Mills theory. 
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normalized as in 



The form of hq for other values of 



Figure 1: Plot of h 

q is similar to this one. The peak moves to the right as q decreases. The locations Zc 
of the peaks are given in Table ^ The plots of hq and hi/q are related by refiection 
around z = 1/2 (see (PUJI ). 
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